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ABSTRACT 

Ideal axisymmetric pulsar magnetosphere is calculated from the standard stationary force- 
free equation but with a new boundary condition at the equator. The new solution predicts 
Ohmic heating. About 50% of the Poynting power is dissipated in the equatorial current layer 
outside the light cylinder, with about 10% dissipated between 1 and 1.5 light cylinder radii. The 
Ohmic heat presumably goes into radiation, pair production, and acceleration of charges - in an 
unknown proportion. 



1. Introduction 

We have shown that ideal pulsars calculated in 
the force-free limit of Strong-Field Electrodynam- 
ics (SFE) dissipate a large fraction of the Poynt- 
ing flux in the singular current layer outside the 
light cylinder (Gruzinov 2011). This result - finite 
damping in an ideal system - is not really that 
unusual. Burgers equation, for instance, with vis- 
cosity +0, dissipates finite energy in infinitely thin 
shocks. 

The standard axisymmetric pulsar magneto- 
sphere features a nearly head- collision of 
Poynting fluxes right outside the light cylinder. 
It is to be expected, although merely by common 
sense, that such a collision should be accompanied 
by damping. 

Here we show that our SFE solution also ob- 
tains from the standard force-free magnetosphere 
equation of Scharlemant & Wagoner (1973), if one 
uses the "correct" boundary condition at the equa- 
torial current layer. 

We propose that the "correct" boundary condi- 
tion at the singular current layer (which now exists 
only outside the light cylinder) is 

B 2 -E 2 =0. (1) 

This condition is Lorentz invariant, comes up in 
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Fig. 1. — Everything is in pulsar units. r s = 0.25. 
Lower Left: Isolines of ^, integer multiples of 
O.lV'o, V'o = 1.44. Upper Left: Poynting flux L 
and the field invariant / = B 2 — E 2 at equator 
vs r. Multiple curves for / are different discrete 
approximations - a rough rendering of numerical 
accuracy. Upper Right: F. Lower Right: Poynting 
flux L through the light cylinder on the field lines 
with z) < tJj. 
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the SFE simulations d, and has a clear physical 
meaning (at equator, the field becomes electric- 
like in order to drive large current). 

We cannot be sure that our proposal works, un- 
til one justifies the full SFE, or just eq.([T]), micro- 
scopically. But conversion of 50% of the Poynting 
flux into the Ohmic power (radiation, electron- 
positron pairs) occurring close to the light cylin- 
der must have consequences for the pulsar phe- 
nomenology, and needs to be studied. 

In §2 we derive the pulsar magnetosphere equa- 
tion and explain how Contopoulos, Kazanas & 
Fendt (1999) solve it. In §3 we put together all the 
equations which are needed to calculate the pulsar 
magnetosphere. In §4 we describe the numerical 
solution and the corresponding physics results. 

2. Ideal pulsar magnetosphere 

Goldreich & Julian (1969) proposed that neu- 
tron star magnetospheres obey the force-free con- 
dition 

/)E+jxB = 0. (2) 

Surprisingly, it turns out that this simple equa- 
tions allows a full calculation of the pulsar mag- 
netosphere (Scharlemant & Wagoner 1973, Con- 
topoulos, Kazanas & Fendt 1999, Gruzinov 2005, 
Spitkovsky 2006). 

For the stationary axisymmetric case, the cal- 
culation is as follows. Using axisymmetry and 
stationarity, in cylindrical coordinates (r, 6, z), we 
represent the fields by the three scalars <f>, if), and 
A, which depend on r and z but not on 8: 

E = -V<^, B = -H 2) 4,i), (3) 
r 

where the subscripts denote the partial deriva- 
tives. 

We plug (3J into ^ and use p = V ■ E and 
j = V x B. We also use the boundary conditions 
at the surface of the star - the continuity of the 
normal component of the magnetic field and the 
tangential component of the electric field. We use 
the pulsar units 

fj, = Q = c=l, (4) 

where \i is the magnetic dipole moment of the star. 
It is assumed that the magnetic field is a pure 

2 See Fig. 5 of Gruzinov (2008). 



dipole near the surface inside the star. The star 
is assumed to be a perfect conductor. f2 is the 
angular velocity of the star. 
We get 

= V, A = A(V), (5) 

where A is an arbitrary function of ij), and we 
also get the "Grad-Shafranov-like" pulsar magne- 
tosphere equation for ip 

(1 - r 2 )AiP - -ip r + FN) = 0. (6) 
r 

Here A = V 2 , and F = A A' , where the prime de- 
notes the ^-derivative. The pulsar magnetosphere 
equation ^ is solved outside the star 

r 2 + z 2 >r 2 sl (7) 

with the boundary condition at the surface of the 
star 

2 

i)= r -, r 2 + z 2 = r 2 . (8) 

The pulsar magnetosphere equation (JSJ con- 
tains F - an arbitrary function of ij), and it is 
not clear how one should solve it. This was ex- 
plained and done by Contopoulos, Kazanas & 
Fendt (1999) (CKF). 

The pulsar magnetosphere equation is elliptical 
both inside and outside the light cylinder, and can 
therefore be solved if some conditions are given at 
all boundaries and if F is known. We first pick a 
trial F. 

Then inside the light cylinder, we have all the 
boundary conditions: (i) we know tj) at the surface 
of the star, (ii) ip = at z = ±oo, (iii) at r = 
1 the boundary condition is given by the pulsar 
magnetosphere equation itself, ip r = F/2. So we 
can find tp inside the light cylinder, say by the 
relaxation method using the variation principle. 

Outside the light cylinder, CKF postulate the 
boundary condition at the equator, 

V>(r,0) = V(l,0), r>l. (9) 

Then, with some boundary conditions at infinity 
(we use %j) = at z = ±oo and ip r = at r — oo), 
and with the same boundary condition at the light 
cylinder, ip r — F/2 at r = 1, one can solve the 
pulsar magnetosphere equation outside the light 
cylinder too. 
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For a generic F, this procedure gives a solution 
with %p(l — 0,z) ^ ip(l + 0,z). But one might 
hope that there is a (unique?) F which gives a 
smooth solution. CKF use a feedback procedure 
— numerical adjustment of F leading to a smooth 
light cylinder crossing. 

In our case, the boundary condition at the equa- 
tor outside the light cylinder is 

(r 2 - 1)(V^) 2 = A 2 , z = ±0, r > 1. (10) 

At the light cylinder, r = 1, this gives 

A(V> )=0, ^ = V(1,0), (11) 

meaning that there is no singular return current 
on the field line tp = ipQ. The only singularity is 
the equatorial current layer. 

3. Pulsar magnetosphere equation 

In summary, we solve the following pulsar mag- 
netosphere equation 

(1 - r 2 )Ai/> - -ip r + AA' = 0, (12) 
r 

V> = 4' r 2 +z 2 = r 2 s , (13) 

(r 2 - 1)(VV>) 2 = A 2 , z = ±0, r > 1. (14) 

The function A{ip) is (uniquely?) determined by 
the continuity at the light cylinder. The Poynting 
power is L = J dip A. 

4. Numerical Solution and Results 

Our numerical procedure is as follows. We pick 
a trial function 

g(r)>0, g'(r)<0, r > 1, <?(1) = 1. (15) 

Instead of the boundary condition (|14j) . we impose 

^(r,0) = V>o<?(r), r>l, ^ = ^(1,0). (16) 

We apply the CKF relaxation method to calculate 
F(ip) (ignoring the singular return current). Then 
we calculate the field invariant at the equator out- 
side the light cylinder 

I{r) = ^(A 2 + (V^f) - (V^) 2 , z = ±0, r > 1. 

(17) 



We then manually chose g(r), so as to make I(r) 
as close to zero as we can at all r > 1. It turns 
out, that 

5(0=0.52+—, (18) 
r 

nullifies I(r) to about the numerical accuracy 0. 

The resulting magnetosphere, the function 
F(tp), the Poynting flux at different radii, and 
the invariant I(r) are shown in Fig.l. 

One gets the spin-down power (the Poynting 
flux at the stellar surface = the Poynting flux 
through the light cylinder): 

u 2 il 4 

£ sd «0.9-^. (19) 

The Ohmic power (the Poynting flux on the field 
lines that cross the equatorial current layer) is 

iohm ~ 0.5L sd . (20) 

The Ohmic power between 1 and 1.5 light cylinder 
radii is 

iohml.5 ~ 0.1i s d. (21) 
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3 This numerical procedure, though formally correct, is 
methodologically inappropriate. In fact, we have 
"cheated". The function g(r) has been read off the SFE 
time-dependent simulation (we then confirm that various 
other profiles g(r) don't nullify I(r) to the accuracy shown 
in Fig.l.). What one really wants is a CKF-type feedback 
loop, which would adjust g(r) so as to nullify I(r). Or, 
maybe, one can enforce the correct boundary condition (114 I t 
throughout the F relaxation. We were unable to develop a 
numerical scheme which would solve the problem (1121 — Hit 
all by itself. 

This 2-column preprint was prepared with the AAS IATgX 
macros v5.2. 
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